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INTRODUCTION 


Concerning  the  stability  of  a uniform  free-free  beam  under 
constant  thrust  with  and  without  directional  control,  the  second 
lowest  branch  of  eigenvalues  turns  out  to  be  the  critical  one  in 
deciding  the  structural  stability  as  it  has  been  reported  by  the 
writer  in  several  previous  papers  [1,2,3].  This  fact  was  not 
mentioned  in  earlier  investigations  [4,5,6].  However,  due  to  a 
numerical  error  as  pointed  out  by  Sundararamaiah  and  Johns  [7],  some 
of  the  writer's  conclusions  also  become  invalid.  This  note  attempts 
to  clarify  some  aspects  of  this  crucial  branch  of  eigenvalues.  First, 
an  analysis  is  presented  showing  that,  other  than  the  first  branch 
of  zero  eigenvalues  (corresponding  to  the  rigid  body  translation  mode), 
there  can  be  only  a set  of  discrete  values  of  the  thrust  which  admit 
zero  eigenvalues.  Stability  data  are  then  presented  which  are  revised 
from  previous  results  with  the  numerical  error  eliminated.  Finally 
conclusions  concerning  the  second  lowest  branch  of  eigenvalues  are 
drawn  from  these  new  evidences. 

ON  DISCRETE  VALUES  OF  THE  THRUST  WHICH  ADMIT  ZERO  EIGENVALUES 

Referring  to  notations  adapted  in  [1],  the  lateral  stability  of  a 
uniform  free-free  beam  under  a constant  thrust  subjected  to  directional 
control  can  be  described  by  the  following  boundary  value  problem: 

D.E.  u"M  + Q(xu' ) 1 + X2  u = 0,  (la) 
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B.C. 


u"(0)  = u-'"  (0)  = 0 


(lb.lc) 


u"(l)  = 0,  u'"  (1)  - K0Qu'(l)  = 0 (ld.le) 

As  mentioned  in  [1],  Q is  the  non-dimensionalized  thrust,  X,  the 
eigenvalue  which  dictates  the  stability  behavior  and  K0,  the 
directional  control  parameter.  A positive  K0  denotes  a rotation  of 
the  thrust  in  the  same  direction  as  the  tangent  at  the  tail-end  and 
a negative  K , in  opposition  to  the  same  tangent. 

First  we  shall  present  here  an  analysis  showing  that,  other  than 
the  rigid  body  translation  mode  for  which  the  eigenvalue  X is  zero 
for  all  values  of  Q,  there  can  only  be  a set  of  discrete  values  of  Q, 

designated  by  Q1 , i=0,l,2 at  which  X take  zero  values  and  that 

Q.'s  are  independent  of  K . For  this  purpose,  set  X=0  in  eqs  (1) 
and  one  has 


D.E. 

u""  + Q (xu ' ) ' = 0 

(2a) 

B.C. 

u"  (0)  = u'"  (0)  = 0 

(2b, 2c) 

u"  (1)  = 0,  u'"  (1)  - 

KQu'(l)  = 0 (2d,2e) 
0 

Eqs  (2)  become 


D.E. 

v1"  + Q(xv) ' = 0 

(4a) 

B.C. 

v'  (0)  = vM  (0)  = 0 

(4b, 4c) 

v'  (1)  = 0,  v"  (1)  - K0Qv(l)  = 0 

(4d,4e) 
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From  eq  (4a)  one  has 


v"  + Qxv  = constant  = c 

The  constant  c must  vanish  due  to  eq  (4c).  In  order  to  satisfy  eq. 
(4e),  one  can  simply  replace  eqs  (4)  by  the  following  system: 

D.E.  v"  + Q[x  - aQe  5(x  - l)]v  = 0 (5a) 

B.C.  v'  (0)  = v'  (1)  = 0 (5b, 5c) 

where 

“e  - 1 + Ke  <«> 

6(x)  is  the  Dirac  delta  function  and  e is  an  infinitesimal  number 
defined  by 

e<5(0)  = 1 (7) 

It  is  an  easy  matter  to  verify  the  equivalence  between  eqs  (4)  and 
(5).  In  particular,  substituting  x = 1 in  eq.  (5a),  one  has 

v"  (1)  + Q(1  - ct0)  v'  (1)  = 0 
which  is  eq  (4e)  by  noting  eq  (6). 

To  solve  eqs  (5),  the  method  of  finite  element-unconstrained 
variations  [2]  is  again  used  here.  It  is  hardly  necessary  to  describe 
the  details  since  they  follow  so  closely  as  in  [2].  Only  the 
unconstrained  variational  statement  is  given  here.  Through  integration- 
by-parts,  it  is  easily  shown  that  the  one  associated  with  eqs  (5)  is 

<$J  = 0 (8a) 

where 

J(v,v*)  = f1  (-v 1 v* 1 + Qvv*)dx  + ea  Qv(l )v*(l ) (8b) 

0 0 

and  v*  is  the  adjoint  variable.  Note  that  v and  v*  are  symmetric  in 
eq  (8b)  indicating  that  eqs  (5)  is  now  a self-adjoint  problem.  The 
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significance  that  e is  infinitesimally  small  (eq  (7))  is  now  apparent. 
As  long  as  a0  is  a finite,  the  solutions  for  Q in  eqs  (8)  are 
independent  of  (or.  Kg  = a^-l).  The  eigenvalue  problem  of  eqs  (8) 
also  indicates  that,  other  than  the  trivial  solution  v=u'=0  (which 
corresponds  to  A=0  for  all  values  of  Q),  the  discrete  eigenvalues  Q.. , 

i =0 ,1,2, are  the  only  solutions.  The  first  eight  of  Q. 's 

obtained  are  listed  in  Table  1. 


TABLE  1.  THE  FIRST  EIGHT  LOWEST  VALUES  OF  THE  THRUST  WHICH  ADMIT  ZERO 
EIGENVALUES  (RIGID  BODY  TRANSLATION  MODE  EXCLUDED) 


i 

0 

1 

2 

3 

4 

5 

6 

7 

Qi/ir2 

0.0000 

2.5677 

9.7218 

21.348 

37.486 

58.156 

83.333 

113.30 

It  might  be  worthwhile  to  mention  a different  way  to  obtain  the 

Q.'s.  Since  the  Q.'s  are  independent  of  K = -1  (afl  = 0)  in  eqs  (5), 

1 1 0 0 

one  then  has*  [8] 

D.E.  v"  + Qxv  = 0 (9a) 

B.C.  v'  (0)  = v'  (1)  = 0 (9b, 9c) 

Eqs  (9)  is  reducible  to  Bessel's  equation  (see,  for  example,  [9])  and 
it  can  be  readily  shown  that  solutions  to  eqs  (9)  are  those  of 

- 1/2 

J2/3(fQ  ) = ° 00) 

*Th1s  is  exactly  the  equations  solved  by  Silverberg  [8],  although  his 
interpretation  of  the  result  was  incorrect  as  it  was  pointed  out 
earlier  by  the  writer  [2]. 
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where  J in  eq  (10)  is  the  Bessel's  function  of  the  first  kind.  The 
first  non-zero  solution  of  eq  (10)  turns  out  to  be  2.5992tt2  which 
agrees  quite  well  with  given  in  Table  1. 

STABILITY  DATA  REVISED 

In  this  section,  the  stability  data  revised  from  previous  results, 
some  of  which  are  misleading  due  to  a numerical  error,  are  presented. 
Figure  1 is  for  the  case  without  control.  The  deformed  parabola 
corresponds  to  the  third  and  fourth  lowest  branch  of  eigenvalues  (for 
the  range  of  thrust  shown).  These  two  branches  now  coalesce  at 
Q = 11.126it2  with  X = 23.10  which  agree  with  the  results  given  in  [7]. 

The  first  branch  is  zero  for  all  values  of  Q and  corresponds  to  the 
rigid  body  translation  mode.  The  second  branch  is  a troublesome  one. 

From  the  present  numerical  results,  for  Q ranging  from  zero  to  IOOtt2, 
say,  the  eigenvalues  of  this  branch  are  too  small*  to  be  distinguishable 
from  the  first  one  in  which  the  eigenvalues  are  supposed  to  be  true  zeroes. 
On  the  other  hand,  the  analysis  of  the  previous  section  has  shown  that 
there  is  no  zero  eigenvalue  solution  other  than  those  of  Q-j . Thus  we  con- 
clude that  this  second  branch  for  Kq  = 0 is  an  "almost  zero"  branch. 

Figures  2 through  5 show  the  first  four  eigenvalue  branches 
(including  the  zero  branch)  for  K0  f 0.  In  the  order  of  the  figure 

*For  larger  values  of  Q (Q  = IOOOtt2,  for  example),  this  branch  appears 
to  be  distinguishable  from  zero  again. 
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numbers  2,3,1 ,4  and  5,  it  is  revealing  to  observe  the  changing  trend 

of  these  curves  as  K varies  from  -1.0  to  -0.05,  to  0,  to  0.05  and  to 

0 

1.0.  The  changing  pattern  of  the  "second"  branch  is  of  particular 

interest.  It  indicates  that,  for  small  |K  | (additional  data  include 

cases  in  which  | | is  as  small  as  0.0001),  the  second  branch  has  a 

divergence  region  in  0<Q<Q1  so  long  as  KQ  is  negative  and  it  becomes 

a stable  region  so  long  as  is  positive.  Hence  the  case  = 0 

appears  to  be  bordering  between  stability  and  divergence  instability. 

In  case  of  divergence  instability,  the  eigenvalue  X indicates  the  rate 

at  which  the  structural  disturbance  will  grow.  Since  the  second 

branch  X is  very  small  for  K = 0,  it  may  be  considered  a stable  mode 

0 

in  the  practical  sense  until  the  flutter  load  Q = 11.126ir^  is  reached. 

However,  unless  Q = 0,  this  small  X does  not  represent  a rigid  body 

rotation  mode  and  thus  must  be  a blending  mode. 

Finally,  Figure  6 shows  the  three  lowest  branches  of  eigenvalues 

for  various  values  of  K and  for  Q between  0 and  4ir2.  It  can  be 

6 

considered  a revision  of  Figure  5 in  [1]. 

CONCLUSIONS 

Concerning  the  stability  problem  of  a uniform  free-free  beam  with 
and  without  directional  control,  we  now  may  summarize  the  following: 

1.  Barring  the  rigid  body  translation  mode,  X can  be  zero  only 

at  a discrete  set  of  Q's  (designated  by  Q-j , i=0,l,2, the  first 

eight  of  them  are  given  in  Table  1).  These  Q-j's  are  independent  of  Kg. 

2.  As  long  as  K is  positive,  there  is  a region  of  oscillatory 
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stability  between  ()  and  Q, . 

0 1 

3.  As  long  as  K0  is  negative,  there  is  a region  of  divergence 
instability  between  and  Q^. 

4.  For  K = 0,  the  second  branch  of  eigenvalues  represents  a 

0 

bordering  state  of  stability  and  instability.  This  branch  may  be 
considered  stable  in  the  practical  sense  due  to  the  smallness  of  the  A's. 
It  is  a bending  mode  nevertheless. 

It  should  be  pointed  out  that  Items  2 and  3 above  are  concluded 
from  numerical  results  only  and  have  not  been  substantiated  by 
analytical  proofs. 

In  retrospect,  the  writer  was  not  justified  to  have  stated  in  his 
previous  papers  that  the  methods  used  by  Feodos'ev  [5]  and  by 
Matsumoto  and  Mote  [6]  had  led  to  incorrect  results.  He  could  only 
state  that  they  did  not  discuss  the  important  second  branch.  However, 
the  Galerkin's  procedure  used  by  Beal  [4],  which  results  in  a second  zero 
branch  for  all  values  of  Q,  Is  still  inconsistent  with  the  analysis 
presented  here. 
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Figure  5.  Four  lowest  branches  of  eigenvalues,  Kfi  = 1.0. 
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ATTN:  DOCU  $ TECH  INFO  BR 

ALABAMA  35809 


12  copies 
CDR 

DEFENSE  DOCU  CEN 
ATTN:  DDC-TCA 

CAMERON  STATION 
ALEXANDRIA,  VA  22314 


.NOTE: 


PLEASE  NOTIFY  CDR,  WATERVLIET  ARSENAL,  ATTN:  SARWV-RT-TP, 
WATERVLIET,  N.Y.  12189,  IF  ANY  CHANGE  IS  REQUIRED  TO  THE  ABOVE. 


